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The Kugo-Ojima color confinement criterion, which is based on the BRST sym- 
metry of QCD is numerically tested by the lattice Landau gauge simulation. We 
first discuss the BRST symmetry on the lattice and the Gribov copy problem. The 
horizon function defined by Zwanziger based on the {/-linear definition of the gauge 
field and that based on log U are compared. The parameter u which is expected 
to be —15? in the continuum theory was found to be —0.7(5^ in the strong cou- 
pling region. There are about 20% difference between the two horizon functions, 
but in the weak coupling region, the expectation value of the horizon function is 
suggested to be zero or negative. 



1 Introduction 

There are essentially two aspects in the manifestation of color confinement in 
QCD. One aspect is the linear potential between quarks, which was conjectured 
by Gribov ^is a consequence of an enhancement of the singularity of the ghost 
propagatoJd due to the restriction on the gauge field A. 

Another aspect is the absence of free single colored particle state in the 
asymptotic Hilbert space, which culminates in the Kugo and Ojima color con- 
finement criterion based on the BRST(Becchi-Rouet-Stora-Tyutin) symmetry: 
i.e. in the Landau gauge, a coefficient in the two-point function produced by 
the ghost, the antighost and the gauge field becomes — 8 h a , where a and b spec- 
ify the color in the adjoint representation. Analytical calculation of this value 
is extremely difficult and so far no verification was performed. 

In 1994, Zwanziger's developed a lattice QCD theory for Gribov ambigu- 
ity. He claimed that, if the restriction to the fundamental modular region is 
achieved, the gluon propagator at momentum should vanish in the continuum 
limit. In the Lattice QCD test, we address the following problems: whether 
the gluon propagator is infrared finite, how singular the ghost propagator is in 
the infrared region, and whether the Kugo-Ojima color confinement criteriorn 
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is satisfied. 



2 The lattice simulation of the Landau gauge QCD and the Gribov 
problem 

2. 1 The path integral formulation of the gauge fixed theory in the presence of 
the Gribov copy (Fujikawa, Hirschfeld) 

First we give a brief review of the path integral formulation of the gauge 
fixed theory by Fujikawa] and HirschfcldS. They discussed a possible situation 
which may give a way out of the gauge fixing degeneracy problem (the Gribov 
problem). Let us define a gauge unfixed partition function as, 

Z= [ dUe-W u \ (1) 



In derivation of the Faddeev-Popov(Faddeev-Popov) formula of the gauge 
f(U) = 0, one considers the determinant function A(U 9 ), 

v ' \ dU 9 dg J V dA9 d 9 J 

where — A„(£7), and one finds that A(U 9 ) is indeed a function of U 9 , and 
looks at the integral on the gauge orbit, U 9 , 

N(U) = [ dgA(U 9 )S(f(U 9 )). 



Obviously N(U) is an orbit function, that is, N(U 9 ) = N(U). At the intersec- 
tion points of the gauge orbit U 9 with the surface f(U) = 0, i.e., g — gi(U), 
Gribov's copies, the above delta function is transcribed to give 

N(U) = J dgY / ^^S(g-g i (U)) = J dgY^sign{A{U 9 ))8{g 9i (U)). 

These delta functions contribute or ±1, and thus, if the orbit function N(U) 
is non vanishing over all orbits, N(U) ^ 0, then in use of the identity 

dgA{U 9 )5{f{U 9 )), 



N{U) 



the standard FP procedure applies, and one factors out the gauge volume and 
obtains the formula for expectation values of functions, F(U), 

/pvm \l JdUA(U)6(f(U))F(U)e-? s W/N(U) 

[f(U))\f gauge- j dUA(U)5(f(U))e-0 S ( U ) /N(U) ' [> 
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This formula allows the BRST formulation as in the following, and also derives 
a natural simulation algorithm of the gauge fixed theory. 



2.2 The simulation algorithm for the gauge fixed theory (Mandula-Ogilvie) 
and the Gribov problem 

Multiplying by the gauge volumes, the denominator and the numerator of the 
formula (Q), respectively, one recovers the Z in the denominator, and has in 
the numerator, 



dU J dgA(U 9 )S(f(U 9 ))F(U 9 )e- f3S{u) /N(U). (3) 
In use of ( |2.lD , one obtains the numerator as 

J dUj2 s *g n (MU 9 *))F(U 9 *)e-P s ( u) /N(U) (4) 

i 

where U 9i is theiJch Gribov copy on the orbit U 9 . This gives us the algorithm 
in the simulation!^ 



ZJ i-, N(U) 

where the last averaging (} is that of simulation, i.e., that with respect to the 
Boltzmann weight e _/3S , and F{U) is a sign weighted average of F(U 9i ) on the 
gauge orbit U 9 . 

However, in case of gauge non invariant function F(U), this formula would 
rather be interpreted as representing Gribov noise. Thus a modification of the 
gauge such that one chooses a unique copy among others on the orbit is favored, 
that is, a new gauge without Gribov copy, and in that case, the above 
formula is useful in practice. 

2.3 The BRST formulation and the Gribov problem 

The obvious standard FP formula allows BRST formulation, 

(F(U))\ f „ = i J d^eMS J cf(U))F(U)e- 0S{U) /N(U), (6) 
Z' = J d»eM5 J cf(U))e-^/N(U), (7) 



where 6 stands for BRST transformation, and the measure is defined as d/i = 
dU dBdcdc, where dB integration is performed on R and dcdc is suitablly 
defined differentiations with respect to Grassmann numbers, c and c. 

ft was shown by Neuberger that if the gauge fixing function f(U) is a 
smooth function of compact variables U, then the expectation value of gauge 
invariant function F(U) becomes an indefinite fornfl, 

(F(U))\ f „ = 2 ( 8 ) 

which implies that all Gribov copies contribute to give total cancellation, in 
other words, the assumption that the above N(U) is non vanishing on all orbits, 
does not hold, and the formula (|J) is totally meaningless. 

An essential point of his argument is as follows. Let us consider a general 
expression 

Z'(t)= ( 'dnexp(tS I 'cf(U))G(U) (9) 



with a gauge invariant function G(U) as the Boltzmann weight. Then one 
finds that 

^r(t) = J d^{5 J cf(U))exp(t5 J cf(U))G(U) (10) 
can be written from nilpotency of 6, i.e., 5 2 = 0, and from 8G(U) = as, 

^(t) = J d^S (( J cf(U))eMtS J cf(U))G(ljfj . (11) 

If W(U, B,c, c) is an analytic function of the compact variables U, then one 
can show that 

J dnSW{U,B,c,c) = 0. (12) 

dZ' 

Thus it follows that —f-(t) = 0. Together with Z'(0) = / • • • dcG(U) = 0, 

one finds that Z'(t) = 0, so Z'(l) = 0. This concludes total cancellation of 
the Gribov copies of gauges given by analytic gauge functions of U. Thus 
one is forced to consider non analytic gauge functions of U as desired gauge 
functions. As a one dimensional U(l) toy example which avoids jl^), one can 
consider U — e lA , A — Imlogt/, where A is not continuous at U = —1, and 
with definition 6 A = c, one finds that 

r /*oo 
dUdBdcdc S (e- B2 c A) = dA dBdcdc c~ B ^ cc ^ 0. (13) 

J J —TV J— OO 

Although this problem is still open in the lattice BRST formulation, use of 
non compact variables in gauge fixing functions may be helpful. 
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2.4 The Landau gauge and the Gribov problem 

Now we focus on the Landau gauge in SU(3) QCD, that is, f(U) — d^A^U), 
where there are some options of definition A^{U) as 

1. J7-linear Ond; A Xtll = -{U x ^ - U\ graceless part, 

2. use of exponential mafl. U x „ = exp^, A^ = -A x ,„ where 
absolute values of all eigenvalues of A Xjtl do not exceed 47r/3. 

In the latter definition, A X ^(U) is not analytic with respect to compact variable 
U Xi/i contrary to the former one. In both cases the Landau gauge, dA 9 = 0, can 
be characterized in use of optimizing functions Fjj(g) of g, such that SFjj(g) = 
for any 5g. 

1. {/-linear definition; F v (g) = ^^L^ J2 X ^ i 1 - W Re trC/ l, M ) • 

2. use of exponential map; F v {g) = \\A 9 1| 2 = J2 x>fl tr (■''•'.■ ,.■>'/ „) , 
where N is the number of colors and d is the dimension. It is noteworthy 
that ||^4 S || 2 is a continuous function of compact variables U 9 in spite of 
non analytic property of A X ^(U 9 ). 

Gribov copy is generic phenomena in both definitions, there exist a lot of 
local minima of Fjj(g) along the gauge orbit U 9 . Thus the naive Landau gauge 
loses its rigid basis both in the theoretical and in the simulation view points 
for examination of gauge non invariant quantities such as gluon propagator, 
ghost propagator, etc. 

The arguments in the preceding subsections formally applies for a cor- 
responding gauge function containing non analytic Heaviside function. The 
gauge-fixing algorithm in the simulation is required to attain the absolute 
minimum of the Fjj(g) along the gauge orbit. But the global minimization 
is difficuJt in general and developing the efficient algorithm is still an open 
problem! 

3 The Kugo-Ojima confinement criterion and the Gribov-Zwanziger's 
theory 

3.1 Kugo-Ojima 's theory 

A sufficient condition of the color confinement given by Kugo and OjimJi is 
that u ab defined by the two-point function of the FP ghost fields, c(x),c(y), 
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and A v {y) 




satisfies u ab (0) = -S ab . 

The corresponding Euclidian expression can be calculated by 



l£yip(*- B ) <tr f X ^D * [-A„A b ]) ) = - ^)u°V). (15) 



where A a is a normalized antihermitian basis of Lie algebra, V a lattice volume, 
and the ghost propagator is given by 



3.2 Zwanziger's theory 

The fundamental modular region A is specified by the absolute minimum 
along the gauge orbits in the Gribov region f2. 



A = {^|||A|| 2 = Min g ||A 9 || 2 }, Acn = {A\-dD>0 , OA = 0} . (17) 



Zwanziger relaxes the periodicity restriction on the gauge transformation 
g, and imposes larger periods than the original. Then some two points in 
the fundamental modular region A may be bridged to be Gribov copies of 
each other, and one of them is not the absolute minimum of the minimizing 
function along the gauge orbit anymore. Surviving points as the absolute 
minimum consist of core region S ( S C A). In the so defined core region 5, 
a horizon function H(U) given below is negative. 

The Horizon function is defined as follows. Let two point tensor G lxvxy ab 




(16) 



be 




xy 



(18) 



Then H{U) is given as 




(19) 



where E(U) reads as follows; 
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1. in JJ-linear version, E(U) = J2i ~k^ e t r ^z> 

2. in A = log U version, E(U) = p— J^tr (A at S*(^/)A a ) , 

I. a 

x/2 



where Ai — adjA t , and S(x) 



th{x/2)' 



Let us define an average tensor G^ vxy be G pLyxv 8 ab = {G^ uxy ), provided color 
symmetry is not broken. One sees that a Fourier transform of the average 
tensor, takes a form 

G,Ap)S ab = Q ^S ab (V - u a \ (20) 



where e = (E(U)) /V., and that it is related with the horizon function as 
(H(U)} 



V 



lim G MM (p) 

p—fQ 



(21) 



He defined the augmented core region ^ = {U : H(U) < 0} n 
( S C $ C f! ). $ and A are qualitatively similar, and he defined the 
partition function in the path integraLin use of the corresponding Lan- 
dau gauge function f-$(U), and concluded^ in the infinite volume limit that 
limy^oo (H(U))/V = 0. Putting Kugo-Ojima parameter as u ab (0) = -S ab c, 
one finds from Q, @, that 

Q + ( d - 1)c -e=(d-l)(c-^)=0 I 

which is called horizon condition. Since we can measure c and e by the lattice 
simulation, we can check to what extent Zwanziger's horizon condition holds in 
our simulation. With respect to the value e/d, note that the classical vacuum 
is characterized by e/d = 1. 

4 The numerical results 

4-1 The formulation of lattice Landau gauge fixing 

We define the gauge field on links as an element of SU (3) Lie algebra as, 

e A ^ = U Xttl , where A^ = -A Xtfl . (22) 
We perform the gauge transformation as e A *^ = g x e A *^ g x +^ and define 
\8A\ = Max x ^ a \dA a x \, \\dA\\ 2 = —l—J^trdAtdA,. 
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The Landau gauge is realized by minimizing ||A 9 || 2 via a gauge transfor- 
mation g^Ug, where g — e £ . In order to obtain e, we switch the following two 
methods, depending on the current value of \dA\ in comparison to some critical 
parameter |9A| cr . 

1. When \dA\ > \dA\ cr , e x = j^j^dA x with suitable parameter 1 < rf < 2.2 

2. When \dA\ < \dA\ cr , e = (—d )t D li (A))~ 1 ridA where 1 < r, < 2 is a 
parameter. 

The restriction to the fundamental modular region is not always achieved. 
But, we observed that the obtained norm \\A\\ is larger or smaller than that 
obtained after the smeared gauge fixingLJ within 1% accuracy. 

4-2 The Kugo-Ojima two-point junction and the ghost propagator 
The FP operator is 

M[U] = -(d-D(A)) = -(D{A)-d), (23) 

where the definition of the gauge field is adopted as U =^,e A , and the lattice 
covariant derivative D^{A) — + Ad(A fJi ) is given in refu. 

The inverse, M^[U] = (M - Mip])- 1 , is calculated perturbatively by 
using the Green function of the Poisson equation Mq 1 = (— d 2 )^ 1 and M\ = 
d^Ad(A^(x)). 

The ghost propagator ( |l~6| ) is infrared divergent and its singularity can be 
parameterized as p~ 2( - 1+a \ where p 2 = ^(4sin 2 -^), {—L/2 < fc M < L/2). 

It depends on (3 slightly, but its finite-size effect is smalEl. These qualitatpe 
features are in agreement with the analysis of the Dyson-Schwinger equation!^. 

We measured the lattice version of |u£(0)| on 8 3 x 16, 8 4 12 4 and 16 4 for 
f3 = 5.5 and 6. When (3 — 6 and the lattice size is small, the Polyakov 
loop distribution deviates from the uniform distribution. In this case, we 
perform the Z% rotation by multiplying the global phase e iz ~ t such that the 
distribution concentrate around one angle, before we measure the Kugo-Ojima 
two-point function. 

We obtained that u^(0) is consistent with — c5£, c — 0.7 in SU(3) quenched 
simulation, (3 = 5.5, on 8 4 , 12 4 and 16 4 . 



8 



1 1.5 2 2.5 



Figure 1: The ghost propagator as func- 
tion of the lattice momentum. The data are 
f3 = 5.5(box) and 6.0(triangle), 16 4 . The 
fitted curve is 1.287/p 2 779 for p = 5.5 and 
1.162/p 2 545 for P = 6.0 (dashed). 



Figure 2: The Kugo-Ojima parameter 
as the function of the spatial extent of the 
lattice aL(fm). The data are (3 = 6 (tri- 
angle), and p = 5.5 (box) 8 4 , 12 4 , 16 4 from 
left to right, respectively. 



4-3 The gluon propagator 

The gluon propagator is infrared finite. We parameterized the zero-temperature 
lattice data using the Stingl's Factorised Denominator Rational Approximant 
(FDRA) method. The effective mass of the gluon in the analysis of 8 3 x 16 is 
found to be about 600MeV. 

The infrared finiteness is in accordance with the Kugo-Ojima color con- 
finement mechanism. As stated in the their inverse Higgs mechanism theorem, 
if we have no mass less vector poles in all channels of the gauge field, A" and 
if the color symmetry is not broken at all, it follows that 1 + u = 0. 

5 Discussion and conclusions 

We performed the first test of the Kugo-Ojima color confinement criterion in 
lattice Landau gauge. We observed that the 8 4 , 12 4 and 16 4 lattice data. The 
data of (3 = 5.5 indicates that u = —0.7, and those of (3 = 6.0 are smaller by 
about 10%. 

In the Zwanziger's theory, the two-point function G ulJ (k 2 ) can be expressed 
in terms .of the Kugo-Ojima two-point function as (|20|). Zwanziger's horizon 
condition!! in the infinite volume limit reads as In terms of the Kugo-Ojima 
parameter c, the left hand side can be written as (e/4) + 3c, and the horizon 
condition is that c = e/d. In the table below e\ and e2 stand for e in our 
16 4 lattice simulation of the first and the second option of the gauge fields, 
respectively. If the gauge fixing is performed so that it brings the configuration 
into the core region or the augmented core region and if the infinite volume 
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limit is considered somehow, then the legitimate check of the horizon condition 
can be done. The core gauge fixing is, however, difficult, and even impossible 
in general, which implies that the core gauge is literarily not the gauge, and 
thus we give the direct results in the table. 






c 


GW,(0) 


ei/d 


e 2 /d 


i/z 3 


5.5 


0.712(18) 


3.14(5) 


0.783 


0.657 


0.78 


6.0 


0.628(70) 


2.88(17) 


0.860(3) 


0.693(1) 


0.86 



Although the renormalization factor of the ghost propagator Z3 cannot 
be fitted precisely, its inverse is numerically close to e\/d. Simulation data 
show in general that when (3 becomes larger, e becomes smaller, while c has 
an opposite tendency. This fact itself does not necessarily disprove the horizon 
condition, but our preliminary data of c which is calculated in the A — log U 
version already gives the zero-intersection of c — e^/d in the increase of (3 from 
5.5 to 6. 

This work was supported by KEK Supercomputer Project (Project No. 99- 
46), and by JSPS, Grant-in-aid for Scientific Research(C) (No.11640251). 
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